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Preface 


This  thesis  is  a continuation  of  previous  work  done  at 
the  Air  Force  Institute  of  Technology  on  Monte  Carlo  tech- 
niques with  the  Weibull  probability  density  function.  It 
is  hoped  that  the  work  of  this  thesis  provides  a basis 
for  more  refined  adaptive  procedures  in  testing  for  Weibull 
pareuneters. 
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GOR/MA/76-1 

Abstract 

An  extensive  Monte  Carlo  effort  is  conducted  to 
sequentially  test  between  two  hypotheses  concerning  the 
scale  parameter  of  the  Weibull  distribution  with  an  unknown 
shape  parameter.  The  methods  are  based  upon  a procedure 
described  by  Cox  in  Sankhya  A,  Vol.  25.  The  first  two 
tests  use  test  statistics  that  are  asymptotically  equivalent 
to  the  likelihood  ratio.  The  third  test  uses  the  likelihood 
ratio  with  the  shape  parameter  estimated  and  with  expanded 
boundaries.  A scale  parameter  of  1.0  was  tested  against 
scale  parameters  of  1.5  and  2.0. 

The  location  parameter  is  zero  in  all  three  tests. 

Error  bounds  were  .20,  .15,  .10,  and  .05.  All  tests  are 
vincensored.  The  three  tests  are  then  run  with  a truncation 
point  of  twice  the  expected  sample  number.  Five  hundred 
Monte  Carlo  repetitions  are  used  for  all  the  above  tests. 
Results  of  the  tests  show  the  actual  alpha  and  beta  errors 
and  the  average  sample  numbers. 

Power  of  the  tr\incated  tests  is  computed  using  two 
hundred  Monte  Carlo  repetitions  for  error  bounds  of  .2 


A MONTE  CARLO  STUDY  OF  COMPOSITE 


SEQUENTIAL  LIKELIHOOD  RATIO 
TESTS  FOR  THE  WEIBULL 
SCALE  PARAMETER 

I . Introduction 


The  Problem 

The  purpose  of  this  thesis  is  to  examine  the  Type  One 
amd  Type  Two  errors  of  sequential  likelihood  tests  for  the 
scale  parameter  of  the  Weibull  density  function  when  the 
shape  parameter  is  unknown.  The  general  method  of  sequen- 
tial testing  using  a likelihood  ratio  test  has  been  sug- 

} 

gested  by  Cox  and  Bartlett.  Cox  and  Bartlett s'  works  are 
an  extension  of  the  sequential  testing  procedure  developed 
by  Wald  (Ref  6:5).  Wald  found  that  the  general  method 
using  a sequential  probability  ratio  test  could  reduce  the 
expected  number  of  required  observations  to  approximately 
one-half  the  size  of  a similar  fixed  sample  test  (Ref  26:1). 

Significance 

The  Weibull  function  has  been  observed  to  be  an 
extremely  flexible  probability  function  because  it  can 
model  components  having  increasing,  decreasing,  or  con- 
stant failure  rates.  This  flexibility  is  achieved  through 
variation  of  three  parameters,  scale,  shape,  and  location. 
These  parameters  will  be  more  fully  explained  later  in  this 


chapter.  Sequential  procedures  have  been  developed  by 
Willicuns  (Ref  28)  to  test  for  the  scale  parameter  when  the 
shape  parameter  is  assumed  to  be  known.  In  cases  where  this 
assumption  is  valid,  the  sequential  probability  test  of 
Willicims  is  quite  useful  in  testing  for  the  scale  paramieter. 
The  difficulty  lies  in  testing  for  the  same  scale  paraimeter 
when  the  shape  parameter  is  not  known.  This  problem  is 
likely  to  occur  in  testing  a relatively  new  item  which  does 
not  yet  have  a well-established  estimate  for  the  shape 
parcimeter.  There  are  methods  other  than  those  proposed 
to  sequentially  test  the  scale  parameter  under  these  cir- 
cumstances. One  possibility  would  be  to  graphically  deter- 
mine an  estimate  for  the  shape  parameter  and  then  use 
Williams'  sequential  probability  test  for  the  scale  param- 
eter. The  reason  for  examining  alternative  testing  proce- 
dures is  to  attempt  to  improve  upon  the  power  and/or 
efficiency  of  existing  tests.  Another  possible  advantage 
is  to  eliminate  the  need  for  cumbersome  graphical  methods 
if  adequate  numerical  methods  can  be  developed.  This  could 
allow  use  of  a computer  routine  to  conduct  the  entire  test. 
Figure  1 shows  some  of  the  flexibility  of  the  Weibull  model 
and  also  the  deinger  of  using  a poor  estimate  of  k in  testing 
for  the  scale  pareuneter. 

Weibull  Distribution 

The  Weibull  distribution  was  proposed  by  W.  Weibull 
in  1939  without  mathematical  foundation.  According  to 
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Soviet  literature,  a "rigorous  mathematical  treatment  was 


done  by  B.  V.  Gnedenko  in  1949  [Ref  7:123]."  For  this 
reason,  the  Soviets  refer  to  the  distribution  as  the  Weibull- 
Gnedenko  distribution.  The  Weibull  distribution  is  con- 
sidered the  most  complex  and  also  the  most  popular  of  the 
commonly  used  paraunetric  family  of  failure  distributions 
(Ref  3:16;  2:46) . 

The  distribution  was  further  publicized  in  1949  as  a 
method  for  testing  the  fatigue  life  of  metal  (Ref  10:14). 

The  distribution  has  a probability  density  function  of 

k9  ^(x-w)^  ^•exp(-  ) x^w  c^O  0,k>O  (1) 

f(x)  = 

0;  elsewhere 

The  three  paraimeters  are  defined  aS: 

6 - the  scale  pareuneter,  which  depends  upon  the 
characteristic  life  of  the  function, 
w - the  location  parameter,  which  tells  the  minimum 

value  of  the  first  failure  time  (guaranteed  life) . 
k - the  shape  parameter,  which  characterizes  the 
failure  rate  of  the  function. 

This  thesis  will  consider  the  location  parameter  to  be  equal 
to  zero.  This  assumption  will  not  effectively  hinder  the 
usefulness  of  the  study  because  the  Weibull  function  is 
primarily  a description  of  lifetime.  Allowing  the  life- 
time to  start  and  possibly  terminate  at  zero  is  reasonable 
for  many  cases.  The  Weibull  distribution  in  this  form, 

(w=0) , is  referred  to  as  the  Two  Parameter  Weibull  function. 
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Another  Insight  into  the  meanings  of  the  k and  6 parameters 
with  reference  to  fatigue  failure,  is  to  consider  k as  a 
function  of  the  mean  ultimate  strength  and  9 as  a function 
of  stress  (Ref  10:32). 

The  mean,  u,  of  the  distribution  equals  QT + 1)  where 
r represents  the  gamma  function.  For  most  values  of  k,  the 
geunma  coefficient  is  nearly  equal  to  one,  therefore  the 
scale  parameter,  6,  can  be  considered  a rough  approximation 
for  the  mean.  The  standard  deviation,  a,  equals 

eir(|  + 1)  - + i)]^^^ 

The  two-pareimeter  cumulative  distribution  function  is: 


F(x)  = 


1-exp [- (x/9) ] 
0;  elsewhere 


x>0  k,e>0 


The  density  function  is 


f(x)  = 


ke"^x^"^*exp[-(x/e)^]  x>0  9,k>0 

0;  elsewhere 


It  should  be  noted  that  the  literature  contains  another 

version  of  the  Weibull  function  in  which  the  scale  parameter 

k • 

is  defined  as  the  value,  -e  . This  changes  the  cximulative 


distribution  function  to 


F(x)  = 


1-exp  t-(xVG)] 
0;  elsewhere 


x>0  G,k>0 


and  the  density  function  to 


I x^“^«expt-xVG] 


x>0  G,k>0 


0;  elsewhere 
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where  G equals  6 . 

Some  confusion  is  possible  as  to  which  scale  parcuneter 
is  being  used.  Test  statistics  employing  both  versions  will 
be  constructed  in  this  thesis.  It  should  be  noted  that  in 
all  cases,  the  scale  parameter  under  consideration  in  the 
test  of  hypotheses  is  0.  Parameter  estimation  will  refer 

A 

primarily  to  0 and  k since  G = 0^. 

To  estimate  k and  0,  the  method  of  maximum  likelihood 
is  used.  This  procedure  maximizes  the  joint  density  func- 
tion or  likelihood  function.  Thoman,  Bain,  and  Antle  have 
provided  an  unbiasing  factor  for  the  maximum  likelihood 
estimate  of  k (Ref  25:11).  Petrick  shows  that,  as  sample 
size  increases,  the  unbiasing  multiplier  approaches  one 
and  is  not  significant  for  tests  with  s^unple  size  as  small 
as  six  (Ref  21:28-30). 

The  methods  used  in  this  thesis  did  not  require  unbias- 
ing of  the  maximum  likelihood  estimates.  As  can  be  seen 
upon  observation  of  the  results,  the  seunple  sizes  were 
usually  large  enough  to  cause  the  Thomas,  Bain,  and  Antle 
unbiasing  factor  to  approach  unity.  No  unbiasing  factor 
was  used  in  the  estimation  of  the  scale  parameter  in  this 
thesis. 

n 

The  expression,  "Z,"  will  represent  Z throughout  the 

i*l 

thesis.  The  maximum  likelihood  expressions  for  0 and  k are 

k0'"^"^  Z x^’'  = 0 (6) 
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(7) 


^ - n In  0 + E In  X.  + e"^  In  9 E x.^ 
d X k 11 

-0^Ex^lnx=O 

where  L(x)  is  the  log  likelihood  function, 

L(x)  »nlnk-nln0  + (k-1)  In  (|)  - E (|)  ’ (8) 

Eqs  (6)  and  (7)  can  be  simultaneously  solved  for  0 and  k 

A ^ 

to  provide  the  maximum  likelihood  estimates,  0 and  k. 

Eq  (6)  can  be  solved  directly  for  0. 


Eq  (9)  can  then  be  substituted  into  Eq  (7)  but  this  does 

A 

not  produce  an  analytic  solution  for  k.  Instead,  Eq  (7) 
must  be  solved  iteratively. 

To  evaluate  the  standard  error  of  the  estimates  of 
0 emd  k,  the  maximum  likelihood  information  matrix,  I, 


3^L(x) 

- E 

3^L(x) 

39^ 

i93k 

3^L(x) 

3^L(x) 

3k30 

- E 

303k 

M — 

(10) 


where  E is  the  expected  value. 

For  the  Weibull  function,  f(x;  8,  k) , this  becomes 
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I « 


k" 

r'  (2.0) 


r'  (2.0) 


1 . r"(2.o) 

T 2 — 


(11) 


(Ref  23:45-46) 


where  T' (2.0)  and  r'*(2.0)  are  the  first  and  second  derive 
tives,  respectively,  of  the  geunma  function.  To  evaluate 
these  two  functions  the  following  procedure  was  used.  The 
diganuna  function,  (i{/)  , is  defined  as 


4i(a)  = r*  (a)/r(a) 

(12) 

iHa)r(a)  = r'(a) 

(13) 

and  by  differentiation 

4)(a)r'  (a)  + i,'  (a)T(a)  = T"  (a) 

(14) 

Substituting  for  F(2)  and  F'(a) 

F(2)  = 1.0 

(15) 

and 

F"  (a)  = i(»^(a)  + tl/’  (a) 

(16) 

From  tables  (Ref  1:268) 

t|)(2)  « .4227843351 

(17) 

!()•  (2)  * .6449340668, 

(18) 

therefore 


r"(2)  = .8236806601  (Ref  5)  (19) 


The  determinant  of  this  information  matrix  becomes  a 
positive  constant.  A, 

A - 1 + r"  (2.0)  - (r*  (2.0)^  » 


1.644934066 


(20) 


The  variance-covariemce  matrix  is  defined  as  1/n  times  the 
inverse  of  the  information  matrix. 

For  the  Weibull  function,  f(x;  G,  k) , the  information 
matrix  is 


I 


^ ^ (In  G+tj<(2)) 

G 

^ (In  G+ij;(2)  A 


where 

A » 1 + i|/’ (2)  + Iln  G + i(/(2)]^ 


(21) 


(22) 


The  psi  function  which  was  previously  mentioned  and  evaluated 
is  also  defined  as 


'l'(a) 


r(a) 


(23) 


and  its  derivative  is 


(a)  = ^ A {In  r(a)] 
da^ 

This  infoinnation  matrix  has  a determinant 


(24) 


A = 4-^  [1  + (2)]  (25) 

which  is  positive  definite  for  all  values  of  G and  k.  The 
log  likelihood  function  for  f(x;  G,  k)  is 

L(x)  = n In  k-n  In  G+(k-l)Z  In  x - ^ E x^  (26) 

G 

(Ref  23:52-54) . 
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Sequential  Tests  of  Hypotheses 

In  general,  a test  of  hypotheses  can  provide  three 
results:  (1)  the  correct  decision,  (2)  a Type  One  error 

which  rejects  the  null  hypothesis  when  it  is  true,  [P(Type 
One  Error)  = a] , and  (3)  a Type  Two  error  which  accepts  the 
null  hypothesis  when  it  is  false,  [P(Type  Two  Error)  = 6]. 
It  is  possible  to  compute  the  probability  of  each  type  of 
error.  It  is  also  possible  to  determine  the  probability 
that  a sequential  test  will  terminate  by  stage  n (Ref  8: 
39-40) . A statistical  hypothesis  is  a statement  about  a 
distribution  of  the  random  variable (s).  If  the  statistical 
hypothesis  completely  specifies  the  distribution,  it  is 
called  simple.  If  the  distribution  is  not  completely 
specified  it  is  called  composite.  .An  example  of  a simple 
hypothesis  is  0=0^  where  the  parameter  0 is  the  only  param- 
eter (Ref  14:268).  The  more  common  example  of  a composite 
hypothesis  is  Q>9^.  This  test  allows  a range  rather  than 
a point  estimation  and,  therefore,  does  not  completely 
specify  the  distribution.  In  this  thesis,  a different 
version  of  the  composite  hypothesis  is  used.  It  states 
0*0j^  where  a second  parameter  (in  this  case,  k)  is  unknown. 

Sequential  Probability  Ratio  Test 

The  Wald  Sequential  Probability  Ratio  Test  (SPRT) 
provides  the  background  for  the  tests  to  be  conducted.  The 
SPRT  primarily  involves  tests  of  simple  hypotheses.  Random 
saunple  units  are  selected  sequentially  and,  after  each 
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additional  observation,  the  test  is  conducted.  Each  test 
results  in  acceptance  of  the  hypothesis,  rejection  of  the 
hypothesis  and  acceptance  of  the  alternate  hypothesis,  or 
no  decision.  If  no  decision' is  made,  an  additional  observa- 
tion is  taken  and  the  test  is  conducted  again.  The  decision 
is  made  by  determination  of  the  test  statistic: 


n 


n f,  (x,)/f„(xj 


i=l 


1'  i"  ‘O'  i' 


(27) 


where  is  the  probability  density  function  assuming 

Hq:  9*9q  to  be  true  and  is  the  probability  density 

function  assuming  the  alternate,  6=6^,  to  be  true. 

The  decision  rules  are: 

1.  If  <_  3/(l-a>,  accept  Hq. 

2.  If  ^ (1-6) /a,  reject  Hq  and  accept  H^^. 

3.  If  6/(l-a)  < (1-6) /a,  take  another  observation. 

a and  6 have  been  previously  defined  under  Sequential 

Tests  of  Hypotheses  and  are  considered  to  be  the  pre-assigned 
risks,  a is  frequently  termed  the  producer's  risk  because 
it  is  the  probability  thdt  the  alternate  hypothesis  is 
accepted  when  the  null  hypothesis  is  true.  6 is  usually 
called  the  consumer's  risk  because  it  represents  the  prob- 
adDility  of  accepting  the  null  hypothesis  when  it  is  false. 

It  should  be  noted  that  the  null  hypothesis  for  the  SPRT 
normally  represents  the  preferable  parameter  value. 

Figure  2 is  a graphical  representation  of  a typical 
untruncated  SPRT  where  a decision  to  reject  Hq  is  made  at 


the  sixth  observation. 
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Tests 

Three  tests  will  be  conducted.  All  have  the  same 
null  and  alternate  hypotheses: 


o‘ 

(28) 

1= 

(29) 

The  values  of  the  par5uneters  Q,  G,  and  k,  will  be  unknown 
and  will  be  estimated  by  their  maximvim  likelihood  estimators. 
An  unknown  parauneter  that  is  not  part  of  the  hypothesis 
(such  as  k)  is  called  a nuisance  parauneter. 

The  tests  will  be  modifications  of  the  likelihood  ratio 
test.  The  reason  for  conducting  a likelihood  ratio  test  is 
that  the  Weibull  function  cannot  be  transformed  into  any 
of  the  standard  underlying  distributions  such  as  the  Normal, 
Chi-square  or  F distributions  which  are  extensively  tabled 
test  statistics.  Adapting  sequential  procedures  to  the 
likelihood  ratio  test  will  provide  an  attempt  to  gain  the 


advantages  described  in  Wald's  SPRT. 


Test  One  will  be  an  application  to  the  Weibull  distri- 
bution of  a sequential  likelihood  test  devised  by  Cox  (Ref 
6:5-12).  This  test  will  use  a test  statistic  based  upon 
9 and  k that  is  asymptotically  equivalent  to  the  likelihood 
ratio  test  statistic.  The  background  and  general  procedures 
for  Cox's  test  will  be  given  in  Chapter  II.  The  acceptance 
and  rejection  boundaries  are  based  upon  Wald's  SPRT  but 
also  include  a coefficient  based  upon  the  information 
matrix.  This  coefficient  will  be  derived  and,  for  Test  One, 
shown  to  be  a constant  larger  than  one.  This  coefficient 
expands  the  boundaries  to  account  for  the  variance-covari- 
ance factors  of  the  parameter  estimations. 

Test  Two  will  be  explained  in  Chapter  IV.  It  will  be 
similar  to  the  first  test  by  also  being  an  asymptotically 
equivalent  test.  The  test  statistic  differs  by  using  the 
general  scale  parameter,  G,  with  the  shape  parameter,  k. 

The  boundary  coefficient  will  be  derived  and  shown  to  be 
a variable.  This  means  that  the  linear  boundaries  of  the 
SPRT  and  Test  One  will  be  replaced  by  non-linear  boundaries. 

Test  Three  uses  the  exact  sequential  likelihood  ratio 
test  statistic  for  G and  k instead  of  the  asymptotically 
equivalent  statistic.  The  variable  boundaries  from  Test 
Two  will  be  used  with  this  test.  Although  Tests  Two  and 
Three  use  the  general  scale  parcimeter  in  the  statistic, 
the  scale  parameter,  0,  is  still  the  hypothesized  test 
parameter.  Test  Three  is  explained  in  Chapter  V. 
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Analysis 

Assumptions.  This  thesis  will  be  concerned  only  with 
Weibull  distributions.  It  is  assumed  that  the  distribution 
to  be  tested  is  known  to  be  Weibull.  No  method  of  testing 
for  the  type  of  distribution  will  be  provided.  The  assump- 
tion setting  the  location  parameter  equal  to  zero  has 
already  been  discussed.  The  scale  parameter,  6,  is  assumed 
to  be  9q  or  0^^.  Data  was  collected  to  show  the  results 
when  this  assumption  is  not  true.  This  analysis  will  be 
given  in  Chapter  VIII.  The  Monte  Carlo  method  of  obtaining 
failure  times  assumes  that  the  simulation  provides  an 
adequate  approximation  of  real  events.  No  samples  will  be 
censored.  Additional  assumptions  concerning  numbers  of  ? 

Scimple  runs  and  the  initial  estimation  of  the  parameter,  k, 
are  embodied  and  specified  within  the  methodology. 

Standards . The  only  criterion  to  be  used  to  judge 
the  feasibility  of  the  tests  is  analysis  of  the  results. 

How  do  the  actual  alpha  and  beta  errors  compare  with  the 
desired  alpha  and  beta  errors?  How  many  samples  must  be 
taken  to  achieve  the  desired  results?  Failure  times  are 
obtained  by  Monte  Carlo  simulation  and  tests  were  conducted 
over  a varied  range  and  combination  of  parameter  values. 

Since  the  true  input  pareimeter  values  are  known  in  the 

simulation,  it  becomes  a matter  of  observing  the  percentage  \ 

of  correct  decisions  and  the  average  S2unple  numbers  required 
to  make  the  decisions.  These  observations  are  then  compared 


14 


with  the  input  alpha  and  beta  errors  and  the  calculated 
expected  S2unple  numbers. 
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II.  Cox  Sequential  Likelihood  Ratio  Test 

This  section  begins  with  a derivation  of  the  Cox  test 
statistic  and  boundary  coefficient  for  continuous  distri- 
butions. A previous  validation  using  a non-Weibull  distri- 
bution will  be  discussed. 

Cox  Asymptotic  Test  Statistic 

The  problem  Cox  was  attempting  to  solve  was  to  formu- 
late a general  sequential  test  of  a parameter  in  the  pres- 
ence of  a nuisance  parameter.  For  ease  of  illustration, 
the  general  parameter  to  be  tested  and  the  nuisance  pareimeter 
will  be  symbolized  by  the  Weibull  parameters  of  interest, 

I 

0 and  k,  respectively.  The  value  of  k is  unknown.  Cox 
replaces  Wald's  statistic  where  k is  known. 


where  is  the  log  likelihood  function^ with  a statistic 

/V 

using  the  maximum  likelihood  estimate  of  k,  k. 


(31) 


This  statistic  is  expanded  about  the  true  point,  (0,k); 


^n  = 


3L^  (x^ , 0 , k) 
30 


, a'^L  (X  ,0,k) 

+ |(9,-0o)  (0,t0Q-20)  


30‘ 

3^L  (x  ,0,k) 


(32) 


! 
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Cox  states  that  the  statistic  is  asymptotically  equivalent 


to  Wald's  statistic  in  Eq  (30)  if,  and  only  if. 


, 3“L  (x  0,k) 
1 n n; 


(33) 


n 36dk 

A ^ 

in  probability.  This  convergence  means  that  0 and  k are 
asymptotically  independent  (Ref  6:5-7).  By  Monte  Carlo 
testing,  Harter  and  Moore  have  shown  that  these  parameters 
appear  to  be  asymptotically  independent  based  upon  sample 
sizes  as  small  as  fifty  with  a fixed  sample  test  (Ref  12:563) 
This  thesis  will  look  at  much  smaller  minimum  Scimple  sizes 
than  the  fixed  sample  sizes  used  by  Harter  and  Moore. 

Using  maximum  likelihood  theory  on  the  terms  of  Eq 
(25),  2 

(x„,0,k) 
n n ~ _T 


30' 


•00 


0 Ln(Xn,0,k) 


nl 


3k‘ 


kk 


0 Ln(Xn,0,k) 


nl 


0k 


(34) 

(35) 

(36) 


where 


log  f(x^,k) 


^00  = 


30' 


(37) 


"kk  = 


3 log  f(x^,k)" 


(38) 


3k' 


"0k  = 


3 log  f (x^;0  ,k) 


303k 


(39) 
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It  should  be  noted  that  the  right  sides  of  Eqs  (37,  (38), 
and  (39)  are  the  definitions  of  the  respective  entries  in 
the  general  information  matrix,  Eq  (10) . 

9 and  k satisfy  asymptotically  the  equations 


1 9 L ( X 6 , 

_ 1 n n;  ' 


k) 


n 


36 


A ^1 

'9k<»-«>  * - 5 9k  ‘ 


(40) 


(41) 


Substituting  these  maximum  likelihood  estimates  into 


Eq  (31)  , the  statistic  to  be  used  becomes 
Z„  - (9j-9„)  nigg  16  - 

The  expected  value  of  n6  is 

E(n0)  = ne 

and  the  variance  is 


(42) 


(43) 


(44) 


vine)  = nljj/llggluj  - ifg-). 

The  stochastic  process  (Z^)  is  a random  walk  with  mean 
increment  per  step  of  9 - variance  per  step 

-1 


of  where 
n 


C = 
n 


1 - 


T^aa  “I 

9k 

99  kk 


(45) 


This  allows  use  of  Wald's  theory  for  normally  distributed 
observations  (Ref  6:7) . 

Cox  Asymptotic  Boundaries 

Boundeuries  are  developed  from  Wald's  original  stopping 
boundaries.  The  Cox  boundaries  become  (Ref  6:8) 
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It  can  be  seen  from  £q  (40)  that  increases  in  the  variance 

coefficient,  C , will  widen  the  distance  between  the 
n 


acceptance  and  rejection  boundaries.  It  is  likewise 

apparent  that  decreasing  the  discrimination  ratio, 

will  have  a similar  effect  by  decreasing  the  absolute  value 

of  the  test  statistic,  Z . 

n 

A Validation  of  Cox's  Test  for 
the  Normal  Distribution 

Cox  attributes  his  test  to  an  earlier  method  by 
Bartlett.  The  main  difference  is  that  Cox  estimates  the 
nuisance  pareuneter  solely  on  information  from  the  random 
variable,  x^,  while  Bartlett  makes  two  estimates  of  the 
nuisance  pareuneter,  kg  and  kj^,  based  upon  the  random  vari- 
able, x^,  and  the  given  value  of  each  of  the  test  parameters, 
and  6^  (Ref  4:239-244).  Joanes  compared  the  general 
methods  of  Bartlett  and  Cox  with  the  more  specific  one- 
sided sequential  T-tests  of  Wald  and  Barnard  for  the  normal 
distribution.  Interestingly,  in  his  initial  publication, 
Joanes  failed  to  use  the  Cox  coefficients  on  the  stopping 
boundaries  and  declared  that  the  Cox  and  Bartlett  tests 
were  not  equivalent.  After  receiving  correspondence  from 
Cox,  a correction  was  published  three  years  later  which 
showed  equivalency  between  the  two  methods.  Joanes  was 
testing  a normal  distribution  with  u=0  and  y*a. 

The  operating  characteristic  function  and  the  average 
sample  numbers  of  the  Bartlett  £uid  Cox  tests  appeared 
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nearly  as  good  as  the  results  of  the  specialized  Wald  amd 
Barnard  tests.  The  primary  difference  between  Cox/Bartlett 
and  Wald/Bartlett  was  for  true  mean  values  near  the  middle 
of  the  range  between  the  hypothesized  values.  In  this 
region  Cox/Bartlett  tests  were  less  powerful  (Ref  15:633- 
637;  16:221) . The  Bartlett  test  was  not  pursued  in  this 
thesis  because  it  is  slightly  more  computationally  involved 
than  the  Cox  test.  Furthermore,  since  it  is  an  equivalent 
test,  similar  results  would  be  expected. 
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III.  Cox's  Asymptotic  Sequential 
Likelihood  Test  with  the  Weibull 
Distribution,  f (x? 6,k) 

This  section  contains  the  mathematical  formulation 
of  the  test,  the  methodology  of  the  procedure  to  be  used, 
and  a discussion  of  its  salient  features. 

Derivation  of  Test  One 

The  general  statistic  to  be  used  was  given  by  Eq  (42) 
in  Chapter  II.  It  is  repeated  here: 

It  can  be  observed  that  the  general  test  statistic  can  be 
used  directly  for  the  Weibull  distribution.  The  value  for 
0 is  obtained  from  Eq  (9)  after  solving  for  k.  Igg  will 
be  defined  in  Eq  (50).  The  two  constants,  6^  and  6^,  will 
be  the  hypothesized  input  values.  It  should  be  noted  that 
contrary  to  the  more  common  practice,  the  alternate  is 
larger  than  the  null  parameter  value.  This  also  reverses 
the  alpha  and  beta  errors  from  their  standard  form  but  is 
consistent  with  the  general  Cox  test  procedure.  The  bound- 
aries are  computed  from  the  general  form,  accept  if 

and  reject  if 

Z > C log  (-r^)  (48) 

n n ’ '1-a' 
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From  Eq  (45)  and  using  maximum  likelihood  estimates. 


C = 
n 


t1 


1 - 


9k 


T ^ T ^ 

•^ee-^kk 


(49) 


with  Igg,  and  Igj^  defined  generally  in  Eqs  (37)  , (38)  , 

and  (39).  These  equations  and  the  likelihood  information 


matrices,  Eqs  (10)  and  (11)  provide  the  solution  for  C^; 


= k^ 

60 


(50) 


igjj  = r'  (2.0) 


= 1 + 


(51) 

(52) 


C = 
n 


1 - 


[r*  (2.0) ]" 

;t7i  , i’"(2.o) 


1 -1 


(53) 


Substituting  the  values  of  r'(2.0)  and  r''(2.0)  obtained 
from  Eqs  (14)  , (15)  , and  (16)  : 


^n  = 


“1 


, (.6449340668) 

" 1.8236806608 


= 1.295466348  (54) 


Methodology 

The  steps  to  be  used  for  a single  test  are: 

1.  Generate  a given  minimum  number  of  failure  times 
using  the  null  value  of  the  scale  parameter  and 
a given  value  of  the  shape  parameter  for  the  two 
parauneter  Weibull  function. 
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2.  Estimate  the  shape  parcuneter  by  solving  Eq  (7) 
using  the  Newton-Raphson  iteration  method. 

3.  Estimate  the  scale  paraunetar  directly  with  Eq  (9) 
using  the  maximum  likelihood  estimate  of  the 
shape  parameter. 

4.  Calculate  the  Cox  test  statistic  using  the  esti- 
mate obtained  in  Steps  2 and  3. 

5.  Multiply  the  constant,  C^,  times  the  standard 
Wald  boundaries. 

6.  Perform  the  comparison  of  the  test  statistic 
with  the  boundaries.  If  a boundary  is  exceeded, 
the  appropriate  hypothesis  is  chosen,  the  number 
of  samples  is  noted,  and  the  single  test  is 
complete.  If  neither  boundary  is  exceeded,  con- 
tinue with  Step  7. 

7.  Generate  one  additional  seunple  failure  time. 

8.  Re-estimate  the  shape  and  scale  pareuneters  as 
in  Steps  2 and  3 using  the  original  plus  the 
added  failure  times. 

9.  Re-accomplish  Steps  4-6  and  continue  with  the 
loop  of  steps  until  a decision  is  made  or  a 
truncation  point  of  four  hundred  is  reached. 

10,  If  the  truncation  point  is  reached,  select  the 

hypothesis  whose  calculated  boundary  is  closer  to 
the  value  of  the  test  statistic. 

The  single  test  covered  above  is  repeated  five  hundred 
times  for  each  combination  of  9q,  9j^,  and  k at  each  input 
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risk  level.  The  percentage  of  times  the  alternate  hypothesis 
is  chosen  becomes  the  alpha  error  for  the  specific  combina- 
tion of  parauneter  values.  The  total  number  of  failure 
s^unples  generated  by  all  tests  resulting  in  acceptance  of 
the  null  hypothesis  is  divided  by  the  nximber  of  acceptances 
to  obtain  the  average  sample  number  to  accept  under 
Similarly  the  average  sample  number  to  reject  under  Hq  is 
computed.  The  number  of  tests  which  were  truncated  are 
recorded  for  acceptance  or  rejection. 

The  steps  and  procedures  listed  above  are  then  repeated 
with  sample  failure  times  generated  with  the  alternate 
scale  pareimeter  in  Step  1.  The  beta  error  is  the  percentage 
of  tests  that  accept  Hq  under  the  alternate  hypothesis. 

The  average ' scunple  numbers  to  accept  and  reject  and  the 
number  of  tests  resulting  in  truncation  are  recorded  by 
the  S3une  procedure  used  under  the  null  hypothesis. 

The  choice  of  four  hundred  as  a truncation  point  was 
felt  to  be  sufficiently  larger  than  the  expected  seunple 
number  in  all  test  cases.  It  would,  therefore,  approximate 
an  untruncated  test. 

Generation  of  Failure  Times 

The  required  number  of  random  variates,  F(Xj^), 
uniformly  distributed  from  zero  to  one,  were  generated 
as  needed  with  the  IMSL  routine  GGUBF  (Ref  15).  The  random 
ntimber  generator  was  seeded  for  each  run  with  a two  digit 


random  number 


These  random  variates  were  changed  into  Weibull  failure 
times,  x^,  by  the  transformation 

» 9‘ t-lnd-F(x^))]^/’'  (55) 

where  F(x^)  denotes  the  cumulative  distribution  function 
of  the  uniform  distribution.  The  entire  Monte  Carlo 
simulation  was  conducted  on  the  CDC  6600  computer  system. 

Newton-Raphson  Procedure 

Given  the  function,  f (k) , for  the  maximum  likelihood 
of  k in  Eq  (7)  and  its  derivative,  f'(k),  set 

f (k)  * k - f (k)/f  (k)  (56) 

Let  Pq  be  an  approximation  to  a solution  of  f(k)=0. 

Generate  the  sequence  [p^]  recursively  by  the  relation 
n=l,2....  Continue  solving  for  p^  until 
(Pj^-p^_j^)  is  less  than  a specified  tolerance  (Ref  20:19-24). 
The  tolerance  used  for  the  computer  program  is  .0000005. 

Eq  (56) , after  appropriate  substitution  becomes 

C 

g(k)  « k + x-(l/n  £ (x  In  x))] 

[nk”^+(k/n) (E  x^  In  x)+E  x^*E  x^  ^ In  x] 

An  initial  estimate  of  k must  be  provided  to  the  Newton- 
Raphson  algorithm.  The  input  value  of  k was  used  for  this 
estimate.  This  value  would,  of  course,  be  unknown  in  an 
actual  test.  It  is  used  in  the  thesis  to  decrease  computer 
time  by  providing  quicker  convergence.  Small  tests  were 
conducted  with  arbitrary  estimates  of  k other  than  the  actual 
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value.  These  tests  all  converged  reasonably  quickly.  An 
upper  limit  of  three  hundred  was  placed  upon  the  number  of 
convergence  iterations  but  the  limit  was  never  reached.  As 
each  additional  observation  is  taken  for  the  sequential 

A 

procedure,  the  value  of  k from  the  preceding  sequential  test 
is  used  as  the  initial  estimate  for  the  next  test.  The 
number  of  iterations  on  the  first  sequential  test  was 
about  ten.  The  subsequent  tests  using  the  previous 
estimate  usually  converged  within  two  iterations. 


t 
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IV.  Cox's  Asymptotic  Sequential  Likelihood 


Test  with  the  Weibull  Distribution, 
f (x;G,k) 

This  chapter  will  explain  only  the  differences  in  the 
derivation  and  methodology  from  that  given  in  Chapter  III. 


Derivation  of  Test  Two 

The  log  likelihood  function  used  as  a basis  for  this 
test  is  given  in  Eq  (26)  . The  general  test  statistic, 

Eq  (42)  , will  be  used  with  parameters  G and  k.  Given  Eq 
(9)  for  the  maximiam  likelihood  of  6 and  the  relationship 

Ai  aV 

G = 9 , 

A 

G = - E x*^.  (58) 

n 

Solving  Eq  (45)  for  C^  by  use  of  the  entries  in  the  infor- 
mation matrix,  Eq  (21)  and  Eqs  (37),  (38)  and  (39): 


(59) 


^Gk  “ " ^ ® '"<2)] 

Gk 


(60) 


li^  - xj  II  + <l^'(2)  + (In  G + i|;(2))2].  (61) 

]c 

Now  using  the  dlgamma  values  given  in  Eqs  (17)  and  (18) : 


[ln(G)+. 422784351]^ 

G'^k^ 


(xt)  fir)  [1.6449340668+(ln(G)+. 422784351)^1 


-1 
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V. 


An  Exact  Likelihood  Ratio 
Test  with  Cox  Boundaries 


This  chapter  will  derive  the  test  statistic  and  list 
differences  from  Test  One. 


Derivation  of  Test  Three 

The  statistic  to  be  used  in  the  sequential  likelihood 
ratio  for  the  G parauneter  with  k replaced  by  k.  Starting 
with  the  f(x;G,k)  version  of  the  log  likelihood  equation, 
Eq  (25), 


Z = Ln(x;Gj^,k)/Ln(x;GQ,ic) 

= [n  In  k-n  In  G.+(k-l)E  In  x.  - ^ E x^] 
1 i 

-[n  In  k-n  In  G_+(k-l)E  In  x.  - J-  E x’^] 
0 ^ 

= n ln(Gj/G^)-(i--  i-)  E x" 


(64) 


This  statistic  will  be  referred  to  as  the  "exact"  test 

statistic.  It  will  use  the  variance  coefficient,  C , 

n 

from  Chapter  IV. 


Methodology 

The  specific  changes  to  the  steps  listed  in  Chapter  III 
are  given  below: 

3.  Estimate  the  scale  parauneter,  G,  directly  from 
Eq  (58). 


1 

! 

f 


i 


1 


29 


VI,  Input  Values  and  Results 

Input  Values 

Input  alpha  and  beta  risks  were  chosen  to  provide  tests 
similar  to  MIL-STD-781B  for  exponential  SPRT's.  For  the 
untruncated  tests,  the  values  of  the  desired  risks  used 
were  .2,  .15,  .1,  and  .05. 

The  values  of  9 were  and  9j^»1.5  and  2.  These 

choices  result  in  the  commonly  used  discrimination  ratios, 

9./0J..  The  SPRT  with  known  k allows  substitution  of  other 

k k 

^1  ^2 

values  of  9 because  of  the  relationship,  9^^  “ ®2 

28:34-35).  This  relationship  does  not  hold  for  Test  One 

but  is  applicable  to  Tests  Two  and  Three.  In  Test  One,  9 

and  k appear  in  forms  other  than  9^^  which  is  the  basis  for 

the  substitution.  In  Tests  Two  and  Three,  the  parameters 

.are  always  of  the  form,  9 , in  both  Z and  C . Normalized 

n n 

time  units  are  assumed  for  9 values. 

The  values  of  the  shape  parameter,  k,  are  .5,  .75, 

1.0,  1.5,  2.0,  and  3.0.  These,  also,  are  representative  of 
frequently  occurring  values  in  actual  Weibull  distributions. 

The  minimum  sample  size  for  the  tests  was  varied 
between  3 and  20.  Since  the  thesis  is  essential  an  initial 
look  at  these  methods,  the  minimum  sample  size  was  cut  off 
whenever  smaller  sizes  appeared  to  be  vinproductive  because 
of  large  output  risks.  Increases  in  the  minimum  sample 
size  were  terminated  when  output  risks  were  close  to  the 
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desired  risks.  The  tabulated  data  will  specify  the  minimum 
sample  size  used  for  each  test  run.  Because  of  the  asymp- 
totic nature  of  the  tests,  minimum  sample  sizes  larger  than 
those  needed  for  an  SPRT  would  be  expected.  Table  I lists 
the  minimum  s^unple  values  used  for  each  test. 

Table  I 

Minimum  Sample  Sizes 

Test  0^ 

1 1.5 

2.0 

2.0 

2 1.5 

2.0 

3 1.5 

2.0 

* Test  run  using  input  beta 
Results 

The  output  tables  for  the  untruncated  tests  are  listed 
in  Appendix  B.  Most  of  the  analysis  can  be  accomplished 
directly  by  reading  the  tables.  This  paragraph  will  note 
general  findings.  The  expected  sample  number,  E(n),  in 
the  tables  is  based  on  6^  and  provides  a basis  of  compari- 
son for  the  average  sample  numbers. 

Error . The  Monte  Carlo  method  assumes  perfect  random 
numbers  and  a perfect  system  model.  This  still  allows 
errors  due  to  the  basic  laws  of  statistics.  These  statistical 
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Minimum  Sample  Size 


5 

3,5 

5* 

5,10 

5,10 

5 

3,5 


= 2*alpha(0=2a) 


1 


errors  are  portrayed  by  the  output  alpha  and  beta  errors. 
These  errors  are  expected  to  decrease  proportionately  with 
l/i^  where  n is  the  number  of  Monte  Carlo  repetitions  (Ref 
24:259).  In  the  case  of  sequential  trials,  this  error 
relationship  is  more  approximate  because  the  number  of 
trials  differs  in  each  run.  The  use  of  five  hundred 
repetitions  did  not  completely  dampen  out  all  the  abbera- 
tions  in  the  alpha  and  beta  errors.  Examination  of  the 
output  tables  in  Appendix  B will  show  occasions  where  the 
output  error  for  a smaller  input  error  is  larger  than  the 
output  error  for  a larger  input  error.  This  is  contrary 
to  expected  results  because  the  smaller  input  risks  have 
wider  boundaries.  An  example  of  this  can  be  seen  in  Test  3, 
Table  III-l,  0j^=1.5,k=.5  and  a minimum  sample  size  of  3.  In 
this  case  the  output  error  for  an  input  of  .15  is  .454  with 
an  average  s£imple  number  of  approximately  25.  The  output 
beta  error  for  an  input  beta  of  .1  is  .504  with  an  average 
sample  number  of  approximately  40.  Most  of  these  variations 
seen  in  the  lower  range  of  the  k values  are  due  to  small 
Monte  Carlo  size.  Central  processing  time,  particularly 
for  low  k values,  was  rather  large.  For  example  with 
0j^«1.5,  k».5  and  input  risks  of  .05,  the  central  processing 
time  for  some  of  the  tests  exceeded  3000  seconds.  It  was, 
therefore,  not  deemed  worthwhile  to  extend  the  trials  to 
obtain  more  consistent  data.  Another  source  of  error  is 
error  within  the  computer.  An  example  would  be  round- 
off error.  For  the  degree  of  accuracy  needed  for  this 
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thesis,  round-off  error  was  not  considered.  Error  due  to 
the  estimation  procedure  will  be  discussed  in  the  next  sub- 
paragraph. 

Estimation  of  Pareuneters.  Since  all  three  tests  used 
the  same  maximiim  likelihood  estimation  procedures,  the 
estimation  results  will  be  discussed  jointly  for  the  three 
tests.  Appendix  E provides  sample  values  of  the  means  of 
the  estimated  parauneters  and  their  associated  average  sauofiple 
numbers.  Estimates  of  k were  consistently  higher  than  the 
real  values.  The  estimates  were  generally  closer  for  0j^=1.5 
than  for  9j^=2.0.  This  does  not  appear  to  be  strictly  a 
function  of  the  average  sample  number.  This  can  be  seen 
by  comparing  estimates  for  different  risk  levels  but  the 
same  6 and  k.  In  these  cases,  a rise  in  sample  size  with 
decreasing  risk  level  does  not  obviously  correspond  to 
improved  estimates.  Estimates  of  k were  slightly  higher 
for  an  input  0=9 as  opposed  to  9=9 q*  The  0 estimates 
were  much  closer  than  the  k estimates.  The  G estimates 
differ  from  their  true  values  because  of  the  bias  of  the 
k estimator. 

A comparison  of  the  expected  sample  number,  E(n), 
with  the  average  sample  numbers  in  Appendix  B shows  the 
average  sample  numbers  to  be  lower  for  Test  One  and  higher 
than  the  calculated  E(n)  for  Tests  Two  and  Three.  Calcula- 
tion procedures  for  E(n)  will  be  set  forth  in  Chapter  VII. 
E(n)  for  these  untruncated  tests  was  based  upon  the  actual 


1 


34 


value  of  k.  The  calculated  E(n)  does  not  anticipate  the 
use  of  minimum  sample  sizes  but  rather  allows  the  test  to 
terminate  after  every  observation.  Minimum  sample  sizes 
are  used  to  increase  the  accuracy  of  the  test  by  d2onpening 
the  affects  of  outlying  data  points  early  in  the  test.  Use 
of  minimum  sample  sizes  will  tend  to  increase  the  average 
seunple  numbers. 

Test  One.  This  test  had  the  lowest  average  sample 
numbers  of  the  three  tests.  For  values  of  E(n)  larger  than 
the  minimum  sample  size,  the  average  sample  numbers  were 
less  than  E(n).  In  the  lower  values  of  k,  the  sum  of  the 
output  risks  was  the  highest.  A bias  to  accept  Hq  can  be 
seen  from  the  data.  Table  B-I-4  in  Appendix  B shows  an 
attempt  to  balance  the  output  risks.  The  input  beta  error 
was  doubled  to  increase  the  output  alpha  error.  Instead  of 
a reduction  in  the  alpha  error,  a decrease  in  the  average 
seunple  number  occurred.  Output  alpha  error  showed  little 
change . 

Test  Two.  This  test  was  strongly  biased  to  accept 
for  all  k values.  The  average  sample  number  to  Reject  Hq 
when  was  true  were  higher  than  the  other  tests. 

Test  Three.  The  least  bias  to  accept  or  reject  was 
seen  in  this  test.  Particularly  good  results  were  noted 
in  the  output  risks  for  low  k values,  6=2  and  a minimum 
sample  size  of  5.  The  output  risks  compared  favorably 
with  the  input  risks. 
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VII.  Truncated  Tests 


The  method  of  determination  of  the  expected  sample 
number,  the  truncation  procedure,  and  the  analysis  of  test 
results  will  be  given. 

Expected  Sample  Numbers 

The  expected  Seunple  number  is  derived  from  the  opera- 
tional characteristic  function: 

“ *hTe).  3h<ei  <“> 

where 


and 


A - (l-0)/a 
B - B/(l-o) 


h(6)  is  the  solution  of 


f (x,0j^) 

_f(x,0Q) 


-|h(0) 


f (x,0) 


1 


(66) 

(67) 


(68) 


P(0q)  » 1-a  (69) 

P(0j^)  - 6 (70) 

Eq  (65)  can  now  be  solved  for  h(0Q)  and  h(0j^)  given  P(0q) 
and  P(0^) . Let 
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(1/0q’')  - (l/Sj^^) 

k InOg/Gj^) 

d 

ln[(l-a)/8] 


hj^  - (In  A)/d 


Now  that  P(0)  Ceui  be  solved,  the  expected  value  of  n,-  E(n), 


E(n|0=[P(0)lnB+(l-P(0) )lnA]/E(Z|0) 

(•»  f(x,0.) 

where  E(Z|0)  = j In  ^ f(x,0)  dx 


= ln(0Q/0^)  + (0/0q)  - fe/0^) 


Eq  (75)  becomes: 


E(n|0) 


- (hg+h^)  P(0) 


S - 0 


Using  Eqs  (69)  and  (70) 


E(n|0.) 


“(ho+hi>  - hg 


S - 0. 


E(n|0j^) 


- 8(hQ+hj^) 


S - 0. 


Eqs  (78)  and  (79)  are  the  expected  sample  numbers  for 
the  SPRT.  For  the  Weibull  distribution  with  0j^>0q,  E(n|0Q) 
provides  a larger  value  than  E(n|0j^).  Since  the  actual  tests 
have  a nuisance  pareuneter  and  are,  therefore,  less  certain 
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than  the  SPRT,  the  larger  valuer  E (n  | 0q),  was  used  to  provide 
a slightly  larger  number  of  observations  before  truncation. 

Truncation  Procedure 

The  method  of  truncation  will  be  to  terminate  the  test 
at  a multiple  of  the  expected  sample  number.  This  thesis 
will  use  a multiple  of  two.  This  method  provides  a straight- 
forward vertical  truncation  that  is  unbiased  by  the  trunca- 
tion. Figure  3 shows  Test  One  truncated  at  2»E(n).  Tests 
Two  and  Three  would  be  the  same  except  the  accept  and 
reject  lines  would  not  have  a constant  horizontal  slope. 


n 


+ 

0 


Reject 

Region 

Continuation 

Region 

Accept 

Region 

2«E(n) 

Fig.  3 Truncated  Sequential  Likelihood 
Ratio  Test 


The  truncated  test  can  be  considered  a compromise  between 
a sequential  test  and  a fixed  sample  test.  Since  k is 
unknown,  it  is  not  possible  to  calculate  the  truncation 
point  prior  to  testing  because  E(n)  is  a function  of  k. 
This  could  be  a disadvantage  in  the  case  of  testing  that 
is  expensive  or  requires  extensive  scheduling.  For  the 
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Monte  Carlo  techniques,  E(n)  will  be  recomputed  with  each 
re-estimation  of  k or  in  other  words,  at  each  observation. 
The  mean  of  E(n)  for  each  of  the  five  hundred  trials  will 
be  listed  for  each  run.  Separate  mean  values  will  be 
listed  for  tests  conducted  under  Hg  and  E(n)g  and 
respectively.  The  minimum  sample  size  was  five  for  all 
tests.  Input  9 and  k values  were  the  same  as  the  untrun- 
cated tests. 


Results 

Appendix  C gives  the  results  for  Tests  One,  Two,  and 
Three.  With  respect  to  the  untruncated  tests  of  Appendix  B, 
results  were  as  expected.  The  average  saunple  numbers  were 
less  and  the  alpha  and  beta  errors  greater.  For  the  larger 
k values,  the  average  sample  nximbers  will  be  greater  than 
2«E(n)  because  the  minimum  sample  size  of  five  was  greater 
than  the  truncation  point.  These'  tests  were  actually  fixed 


sample  tests. 


Table  II  lists  the  other  changes  that  were  noted  by 


truncating . 


Table  II 


Conipagison  of  Truncabed  and  Un truncated  Tests 


Test 

M 

<D 

Remarks 

I 

1.5 

No  change. 

2.0 

Improved  output  risks. 

II 

1.5 

No  change. 

2.0 

Less  bias  to  accept  Hq,  improved 

g risk. 

III 

1.5 

No  change. 

2.0 

Caused  a bias  to  reject. 
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VIII.  Power 


The  power  of  the  test  is  the  probability  that  0^ 
will  be  rejected  when  0 is  the  true  parauneter  value.  When 
0=0Q,  the  power  equals  the  alpha  error.  When  0=0^^,  the 
power  equals  1-S.  For  the  SPRT,  Wald  has  devised  a 
computational  method  of  determining  power.  Eq  (65)  is 
the  power  equation.  A more  complete  description  can  be 
found  in  almost  all  discussions  of  the  SPRT.  A recommended 
description  is  Mood  and  Graybills’ (Ref  19:388-391) . Because 
of  the  differences  of  the  sequential  likelihood  ratio  test, 
the  computational  method  was  not  used.  Instead,  a Monte 
Carlo  test  was  conducted  with  fourteen  input  values  ranging 
from  .1  to  2.6  for  the  discrimination  ratios  of  1.5  auid  2.0. 
The  same  six  values  of  k that  were  used  in  .the  prior  tests 
vere  repeated.  The  tests  were  truncated  at  2*E(n).  The 
number  of  trial  runs  for  each  combination  of  inputs  was  two 
hundred.  Input  risks  were  limited  to  .1  and  .2.  The  power 
was  determined  by  svibtracting  the  output  beta  error  from 
one.  Figures  4 and  5 provide  graphical  representations  of 
the  power  of  the  three  tests  for  the  specified  sets  of 
inputs.  These  figures  should  not  be  construed  to  represent 
the  power  relationships  between  the  tests  for  all  inputs. 

The  t5djulated  power  for  all  the  tests  appears  in  Appendix 
D.  Generally,  all  three  tests  showed  typical  power  curves 
except  that  the  errors  at  9^  and  0^  were  greater  than 

designed  values  at  low  k values. 
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IX.  Conclusions  and  Recommendations 


Conclusions 

It  cannot  be  summarily  stated  that  any  one  of  the 
three  tests  is  superior  to  the  others  in  all  facets.  It 
can  be  said  that  all  three  tests  showed  similar  results  to 
those  of  the  SPRT.  Test  One  was  notable  for  having  the 
smallest  average  seimple  numbers.  Test  Three  showed  the 
least  bias  in  the  output  risks.  These  two  tests  are 
probably  the  best  candidates  for  further  development. 

These  basic  tests  provide  a basis  for  trying  adaptive  proce- 
dures to  improve  the  results.  Any  of  the  three  tests  could 
be  used  in  their  present  form  for  the  higher  k values. 

Recommendations 

• Adaptive  procedures  are  quite  extensive  in  sequential 
testing  and  could  provide  significance  improvements  to  the 
tests  in  this  thesis.  One  adaptive  procedure  by  Harter 
and  Moore  that  is  similar  to  Test  Three  has  been  published 
(Ref  11:100-104) . This  test  used  Wald  boundaries  with  k 
corrected  for  bias  and  had  smaller  average  sample  numbers 
than  Test  Three.  Other  adaptive  procedures  that  could  be 
tried  include  sequentially  decreasing  the  boundaries, 
decreasing  the  boundaries  if  the  k estimate  are  greater 
than  1.0  or  1.5,  experimenting  with  different  combinations 
of  input  risks,  and  developing  other  functions  of  0 to 


use  as  the  asymptotic  test  statistic.  Ghosh  states  that 
the  given  approximation  for  the  boundary  coefficient, 
is  essentially  a first  approximation  and  that  more  refined 
approximations  may  be  possible  (Ref  8:333).  The  two  most 
prominent  objectives  would  be  to  improve  the  output  risks 
for  the  low  k values  and  to  increase  the  output  errors  up 
to  approximately  the  input  errors  for  the  higher  k values 
in  order  to  lower  the  average  sample  numbers.  The  first  of 
these  objectives  is  obviously  the  more  difficult.  Improving 
the  k estimation  procedure  might  provide  improvement. 

Another  estimation  procedure  which  could  be  tried  has  been 
published  by  Wingo  (Ref  29;  Ref  30) . This  previously  men- 
tioned unbiasing  factor  of  k could  also  be  tried  with  these 
tests.  It  may  be  possible  to  combine  tests,  using  one 
test  statistic  and  boundary  to  accept  and  another  to  reject. 
These  recommendations  do  not  exhaust  the  list  of  possible 
adaptive  procedures.  They  are  offered  only  in  the  hope  of 
inspiring  further  invesigation. 
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APPENDIX  A 
COMPUTER  PROGRAM 


C SEOUENTIAL  LIKELIHOOD  RATIO  TISTT  WITH  OPTION  TD  TRUNCATE 
C PROGRAM  RirJIRES  ATTACHMENT  OF  IHSL  LIBRARY  TO  OBTAIN  RANDOM 
C GENEPaTING  FUNCTION  GGU9F 

C DISTINGUISHES  BtTHcEM  HtIBULL  DISTRIBUTIONS  WITH  RESPECT  TO 
C SCALE  PARAMtTt®.  UNKNOWN  SHAPE  »ARAMETuR  REPLACED  BY  ITS 
C MAXIMUM  LIYELIHOOP  ESTIMATE*  MONTE  CARLO  SIMULATION 
C TEST  l-ASY**PTOTIC  SLRT-COX  ( IPS  T)  METHOD  WITH  THETA  STATISTIC 
C TEST  2-ASVHPTOTIC  SLRT  USING  COX  METHOD  WITH  G STATISTIC 
C TEST  S-EXACT  SLRT  WITH  G STATISTIC  USING  COX  30RY  COEFFICIENT 
C ••  DEFINITIONS  •* 

C ISEEO-SriO  FOR  RANDOM  NUMBER  GENERATOR 

C lA  - PARAMETER  TO  CHANGE  INPUT  AL®HA  AND  SETA  ERRORS 

C IV  - PARAMETER  CHANGING  INPUT  THETA  TO  ALTERNATE  HYPOTHESIS 

C IP-oARAMETER  TO  CHANGE  SHAPE 

C ASM  - AVERAGE  SAMPLE  NUMBER 

C (A)  - ACCEPT 

C (R)  - REJECT 

C ALPHA(IA)/BETA-IMPUT  RISKS 

C RALPHA/P3ETA  - COMOUTEO  OUTPUT  pRPOR  ( AL®HA/RE’' A) 

C NTRUNC  - TRUNCATION  POINT  FOR  SEQUENTIAL  OBSER'/ATIONS 
C MULT-TRUNCA7IOM  MJLTIPLE  OF  EIN) 

C NROBSA/NROBSF-  NUMBER  OF  03S  UNTIL  ACCE®r /REJECT 
C NRACPT/NRRJCT-NJMBER  OF  TE^TS  RESULTING  IN  ACCE’T/RE JECT 
C NT RU NA C /NT PU NR J -MR  OF  TESTS  RESULTING  IN  TRUNCATION 
C IQ  - CUMULATIV*^  DBS  TOTAL  FOR  EACH  RUN 
C TOL-  tolerance  pOR  NEWTON  RAPHSON  ITERATION 

C SMALL/RLARGF-TOLERANCES  TO  ELIMINATE  DATA  THAT  MIL-  ABORT  RUN 
C NUMORO  NR  OF  BAD  DBS  EXCEEDING  SMALL/RLftRGE  LIMITS 
C 101  - MINIHiJN  NR  OF  OBS  FOR  EACH  RUN 
C NSAMP-NR  OF  rums  PER  TEST 
C •* 

PROGRAM  HOFFKINPUT, OUTPUT! 

OIMEMSION  X(40  3! ,FX <40 0!  , N TRUNAC (4 , 2) ,NTRUNRJ <4» 2 ) 
niMENSION  ASNA(4,2) ,ASNR<4,2) ,NUMDB3{4,2) 

DIMENSION  NR09SA  (4, 2) ,NROBSR (4,2) 
riHEMSION  RNR0BSA(4,2) ,RNROBSR<4,2) 

DIMENSION  PALPHA<4!  ,R3ETA<4) 

DIMENSION  ALPHA (4) ,RKI(6) 

DIMENSION  RKK(5t’:),EcTHETA  (Fu  0)  , EKHE  AN  (4 ,2)‘,  THETAM  (4,  2) 
DIMENSION  RNTRNC(500) 

niMENSIO!'  VARK(4,2)  ♦VART(4,2!  ,VARN  ( 4 ,2)  ,PNTRMN<4 , 2) 

DATA  ALPHA (1)/. 20/, ALPHA (2! /, 15/,AL®HA < T)/. 13/, 

DATA  NSAMP/S0C/,THETA0/1.0/,NTP.UNC/40  0/,AL°HA('4)/*3  5/ 

DATA  TOL/.0a00a05/,SMALL/l.C:-050/,RLARGE/2300./ 

DATA  RKKl)  /.5/,RKI  (2)/.75/,RKI  (3)  / 1. 0/, R<I  ( 4) /I. 5/ 

DATA  RKK5) /2*(3/,RKI<6)/3.  C/ 

DATA  THETA1/1.3/,IQ1/05/,ISE£0/19/ 
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DATA  IAl/l/,Ifl2/4/,IA3/l/,IPl/l/,IP2/5/,I?3/l/ 

DATA  MULT/0  2/ 

C 00  LOOP  FO?  CHANGING  SHA^E  ^ARAMETE^  CO 

DO  111  IP*IP1,IP2,IP3 
RKO=R<I(IP» 

RKOINrf=l./RKO 

C LOOP  FOR  CHANGING  ALPHA/nETA  LEV'L 
DO  1 IA=IA1,IA2,IA2 
SETA=ALPHA(IA) 

AsALOGCd.-HETA)  /ALPHAtIA)) 

BaALOG(B£TA/Cl,-ALPHA(IA) » ) 

C CALCULATE  NUMERATOR  OF  E(NI  IF  TRUNCATING  ON  '1JlTI®L£  OF  E(N) 
RNUMEsALPHA(IA) ♦(AL0G((1-ALPHA(IA) ) /BETA) f A) -A.06 (i-ALPHA 
C(IA))/3ETA 

C LOOP  FOR  CHANGING  INPUT  THETA  TO  ALTERNATE 
no  2 IV=1,2 

IFCIV.E0,1)THE7A=THETAQ 

IF(IV.F0.2)THETA=THETAl 

NTRUNACCIA,IV)=0 

KTRUNRJCIAtIV)*: 

NRAC«»T  = 0 
WRRJCT=0 
NUMOBOf IA,IV)=? 

N?.ORSA(IA,IV)=0 
NR033R<IA,I'/)=0 
C MAKE  NSAMP  Rl.'NS  . 

DO  3 IN=1,NSAMP 
23  1*'*I01 

C LOOP  FOR  GEMFRATING  INITIAL  131  RANDOM  VARIATES  ’ 

00  4 10=1,101 

FX«1D)=GGUPF<ISEED) 

04  CONTINUE 

0 LOOP  TOCOMPUTE  X-VALUES  (FAILURE.  TIMES)  FROM  FX  VA.JES 
00  5 IE=l,ni 

X(IE)»THFTA*«C-AL0G(1.-FX(IE) ))**RKCINV) 

05  CONTINUE 

C IF  LARGEST  X=0,  THROW  OUT  SArtPLE  AND  TAKE  A NEW  ONE 
SUMH>Q.OCOOO 
DO  6 iz=i,n 
SUHM*X(IZ)  ♦•SUMM 

06  CONTINUE 
IF(SUMM.GT.3.)G0  TO  40 
NUHOBOCIA,IV»*NUMOROaA,IVM-l 
GO  TO  99 

C ITEPATIVE  NEWTON-RAPHSON  METHOD  TO  CALCULATE 
C MAX  LIKELIHOOD  ESTIMATE  OF  Kt 

C INITIAL  VALUF  OF  K1  IS  KO 

40  RKl*R<a 

C LOOP  FOP  PERFORMING  ITERATIONS 

41  DO  7 IF=1,300 
SUN 1=3. 

SUM 2=0. 
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SUM3=0. 

SUH4*0 . 

C LOOP  FOR  CALCJLAtlNG  SUMMATIONS 
DO  8 IG=1,I0 
C SUM  OF  X(I)**<1 

T1=X{IG)**R<1 

SUM1*SUH1*T1 

C SUM  OF  LM(X(D)  ^ 

T2sAL3G<X(IG)> 

SUM2sSUM2fT» 

C SUM  OF  (X(n**Kl)*LN(X(D) 

C SUM  OF  (X<I) »*<1)»LM(X(I)) 

SUM3»SUM3*T1*T2 

C SUM  OF  (X(T)**<1)»((LN(X(I)))**2 
FUH4=S'JHt>  + Tl*(T2**2.) 

08  CONTIMUF 

C CALCULATE  NUMERATOR  OF  F(X1)  AND  PART  3F  DENOMINATOR 
Ml»IO*SUHi  - 

W2sn*S'JM 

C CALCULATE  OEL'OMINATOR  OF  F(K1I 
H3*H2-SUf"l*S'JM2 

C IF  OENOMINATCR  = C,  THROW  OUT  SAMPLE  ,TAKE  NEW  ONE, CONTINUE 
IF  CA8S(W3)  .LT.SMALDGO  TO  1 
M4*H3»W3 

C CALCULAT*;  pIPST  DERIVATIVE  OF  F(ri) 

FBK=l.-(  (W3»W2-(Wl*CIQ*SUM4-3UM3*SUM2n  » /(W4)  ) 

C IF  FPK  s a , 00  SAME  AS  WHEN  W3=0* 

IFCA33(FP!<)  .LT.SMALDGO  TO  351 
C CALCULATE  F(Ki) 

FKaR<l-(Wl/W3> 

C CALCULATE  GK  FDR  MFWTON  RAPHSON  MSTHOO 
GKsRKl-(FK/FP<l 

C IF  MEM  VAL’JF  OF  K1  IS  TOO  LARGE,  THROW  OUT  SAMP.E 
IFtGK.GT.RLARG^lGO  TO  461 

C IS  NEW  VALUE  LESS  THAN  TOL  ? IF  YES, DEPART  LOOP,  OTHERWISE 
C CONTINUE  WITH  NEW  <1  » GK 
T£ST*ARS(GK-RK1) 

TFCTEST.LE.TODGO  TO  63 
PKlsCK 

07  CONTINUE 
GO  TO  64 

C LOOP  FOR  CALCULATING  SUM  OF  XCI)**K1 

63  R«ClsG< 

64  SUMOCxQ.CQCiI 

00  9 IHsl,n 

SUM03  = SUM0  0fXaH)**RKl 

09  CONTINUE 

C ESTIMATE  THFTA  (ETHETA) 

C IF  USING  TEST  ONE  USE  FOLLOWING  THETA  ESTIMATE 
ETH£TA=  (SU»«Ofl/IO)  ••  (1/RXl) 

C IF  USING  TPST  2 OR  3 THE  FOLLOWING  ETHETA  IS  THE  ESTIMATE  OF  S 
ETHETA=SUMaa/ID 
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C CALCULATE  T-TST  STATISTIC  CZ) 

C IF  USING  TEST  ONE  USE  FOLLOWING  STATISTI3 

Z=CTHETA1-THETA0)*(R<1**2)  *10* (ETUI T A-, 5* (THETA3 4-THETA  11) 

C IF  USING  TEST  TWO  USE  FOLLOWING  STATISTI3 

Z=(  nH£Tfll**R<l)-(THETA0**RXl))*I3»  (£THETA-,5*(  (THETA0**R<1) 
C4-(THETA1»*RK1)))  /<(ETHETA)  •*2.) 

C IF  USING  TEST  THREE  USE  FOLLOWING  STATISTIC 

ZsIQ^RKl’ALOGCTHETAO/THETAl) K ( l./THETAO »»R<1) - { l./THETAl 
C**R<1))*SUN1 

C CALCULATE  OOUNOARr  COEFFICIENT  (C) 

C IF  USING  TEST  ONE  USE  FOLLOWING  COEFFICIENT 
C=1.29546-6T48 

C IF  USING  TEST  TWO  OR  THREE  USE  FOLLOWING  COEFICIEMT 
TERH1=IAL0G(ETHFTA)  )**2•.^.8U556867•AL0G(ETHETA^ 

TERM2=l,-(  (TERN1+.17B7V6594)/  <TZRNH-l.B23*i83  3oOS)  ) 

CS1./TERH2 

C CALCULATE  OOUNOARIES  (ACPT/RJCT) 

AC«»T=C*3 

RJCT»C*A 

C CALCULATE  TRUNCATION  POINT  IF  TRUNCATING  ON  A 1ULTIPLE  OF  E{N> 
NTRUN;=RNUHE/(RKl*ALOGCTHETAl/rH£TAO)-( (THET Al/ THETAO) **R<1 
04-1)4-1 

NTRUNC=HULT*NTRUNC 

C CONDUCT  TEST 

IFCZ.LT.7CPT)G0  TO  81 
IFfZ.3T.RJCT)G0  TO  82 
IFI  CIO)  .Gc.NTRUNOGO  TO  83 

80  10=104-1 

C OBTAIN  ONE  HOPE  FAILURE  TI^E  AND  ITERATE 
FX(IO)=GGUeF(ISE£D) 

X(I(3)=THETA*((-AL0G(1.-FX(IQ)  ))  **R<uINV) 

GO  TO  41 

C INCREHENT  ACCEPTANCE  COUNTER 

81  NRACP7  = NRACPT4.1 
NR08SA(IA,IV)  = NR09SA(IA,Itf)  HO 
GO  TO  303 

C INCRFMFNT  REJECTION  COUNTER 

82  NRRJCT  = NPRJCT4-l 
NR09SR(IA,IV)=NROBSR(IA,IV)4-I0 
GO  TO  3C3 

C TRUNCATION  Tf  ST 

83  TRUNCV  = < ACPT4-RJCT)/2. 

IF(Z.LT.TRUNCV)GO  TO  84 
NTRUNRJ(IA,IV)=NTRUNRJ(IA,  IV)  4-1 
60  TO  82 

84  NTRUNAC(IA,IV)=NTRUNAC(IA,IV)  4-1 
GO  TO  81 

C TAKE  ANOTHER  SAHPLE 

C COUNT  BAD  OBS  AND  REPLACE  SAH»LE 
GO  TO  61 
GO  TO  61 

61  NUH0B0(IA,IV)  = NUN090(IA,IV)4-1 


L. 
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99  IFtNUinaoCIA,!'/)  .GE.30)  GO  TO  101 
GO  TO  23 

C CALCIILATE  HEAM  AND  VARIANCE  OF  K, THETA, ANO  ASM  ESTIMATES 
3C3  RKKCIM)=PK1 

PNTRMC ( I W) =NTRUNC/ MULT 

C IF  USING  TEST  1 USE  COLLOMING  TO  FINO  MEAN  OF  THETA  ESTIMATE 
EETHETA(IH>=£THpTA 

C IF  USING  TEST  2 OR  3 USE  FOLLOWING  TO  SET  MEAN  THETA  ESTIMATE 
fcFTHETA(IU) sETHETA**(1/RK1) 

rnsamr=nsam» 

D3  CONTINUE 
?UHN=0.C 
SUMT*0.0 
SUHK=0.0 

OO  13  IMWaltNSAMP 
SUMN«SUMNfRNTRNCawW) 

SUMK=RKK(IMH)fSU*<K 
SUMT*£ETHETA  tlWW ) fSU*<T 

13  COMTIMUE 

RNTPM'4(IA,IV)=SUMN/NSAMP 
EKMEAN  CIA,TV)sSUHK/NSAMP 
THETAM(IA,IV)sSUMT/NSAMP 
SUMVNsO. 

SUMVK»0. 

S»JMVT=0. 

OO  14  IM=1,NSAMP 

RNTRNC(I**»  = (RNTRNC(IMI-RNTRMN(IA,IV))  **2 
RKK(IM)  = (R«(IH)-FKHfAN(IA  ,I7))**2 
EETHETA(IM)  = (EETHETA(Irtl-THETAMaA,  If)) 

SUMVN=3UMVN»>9NTRNC  ( IM) 

SUMV<=SU“!V'<i>RK<(IM) 

SUHVTsSU^'VT+EETHETA  (IM) 

14  CONTINUE 

VARN(IA,TV)sSUMVN/(NSAMP-l.) 

VARK(IA,IV)sSUMy</(NSAMP-l.) 

VART(IA,IV)*SUMVT/(NSAMP-1.) 

FINO  ASN,  ALPHA,  AND  PETA  ERRORS  FOR  EACH  RUN 
CONVERT  INTEGERS  TO  REAL  NUMBERS 
PNRAC’TaNRACPT 
P.NRRJ3T=NRRJCT 

RNR  OPS  A < I A , I V)  = NRO  9S  A ( I A , I V) 

RNROaSRCIA,IV)  =NROBSR (I A, I V) 

0BSAVSA=FNR09SA (IA,IV)/RNRAC’T 
IFfIV.£0.2)G0  TO  91 
RALPHA (IA)=RNRRJCT/RNSAMP 
91  ASNA(IA,IV)sOaSAVGA 

orsav3R*rnrobsr(:a,iv)/pnrrjct 
IFdV.tO.DGO  TO  92 
RBETA ( lA )sRNRACPT/RMSAMP 
.92  ASNR(IA,IV)3QBSAVGR 
32  CONTINUE 
01  CONTINUE 
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C P^INT  RESULTS 

oRINT»,“  SHAPE  PARAMET£R(<)  = ",RKO 

C OUT»»UT  DATA 

95  PRINT  96 

96  FORMAT  (IH  , TV, “INPUT  0UT»JT  »*»*45'^**** 

CTRUNCATIOM  ESTIMATES") 

PRINT  9'? 

97  FORMAT (IH  »2Y, “THETA  KC  ALPHA  AL»HA  BETA  HO  HI 

C HO  HI  KO  THETA3  K1  THE 

CTA1“) 

PRINT  3B 

98  FORMAT (IH  HO  HI  J3ETA  ACPI  RJCR  ACPT  RJOT 

CACC  REJ  ACC  REJ  MEAN  VAR  MCAN  VAR  MEAN  VAR  MEAN 
CVAR“) 

00  12  IA  = IA1,IA2,IA-! 

PRINT  3?,THETA3,THETA1,R<0, ALPHA(IA) ,RALPHA(IA» ,R3ETA(IA), 
8ASNA(I  A,l)  , ASNRdA,  1)  ,ASNA  (lA  ,2)  , ASNR  (lA , 2)  , NTRJNAC  (I A , 1)  , 
9NTRUNRJ(IA,1)  ,NTRUNAC(IA,2)  ,NTRUN.RJ(IA,2)  ,E<MEAN(IA,1)  ,7AR< 
C(IAtl)  ,THETAM<IA,1)  ,VAP.T  ( I A ,i ) , EKME  AN  ( I A , 2 ) , V ARK  ( I A ,2) 
C,THETA1(IA,2) ,VART(IA,2) 

85  FORMATdH  , ^3. 1 1 IX,  FT.  1»  IK  , F<» . 2, 2X , r 3. 2, 2X  , . 3 , 2 X , ^4.  3, 1< , 
CF5.1,1X,F5.1,1X,P5.  1,1X,F3.1,1X,I3, 1X,I3,1X,I3,1X,I3,3X,F5. 
C2,2X,F5. 2,2X,F5.2,2X,F4,1, 2X, F5 . 2, 2X, F4. 1 , 2X,  F5 . 2 , 2X,F6. 2) 

12  CONTINUE 

PRINT*, “ TRUNCATION  POINT=  “, MULT,“  *E(N)“ 

PRINT  20 C 

200  FORMAT (1H,2X, "MEAN/VAR  OF  ECN)  HO  HI 

C ALPHA") 

00  113  IAsIAl,IA2,IA3 

PRINT  52,RNTRMN(IA,1) ,VARN (IA,1) ,RNTRMN{I A, 2) ,VARNCIA,2) , 
CALPHA(IA) 

62  FORMAT (1H,25X,F?.1,2X,F7.2,4X,F5.1,2X,F5.2,14X,F3.2)‘ 

113  CONTINUE 


111  CONTINUE 
101  STOP 
END 
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Appendix  B 

Output  Risk  and  Average  Sample  Number 
for  Truncation  Point  of  400 

Risks  - desired  alpha  and  beta  errors. 
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Numbers  in  parentheses  behind  average  sample  numbers  are 
the  number  of  tests  resulting  in  truncation. 

E(n)  - expected  sample  number  using  input  parameter  values. 

Tables  are  based  upon  500  runs. 

Table  B-II-4  uses  an  input  beta  risk  equal  to  twice  the 
alpha  risk. 
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Output  Risks  and  Average  Sample  Number,  Test  One 
, 0,=1.5  Minimum  Sample  Size=5  Asymptotic  0 Statistic 
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Appendix  C 

Output  Risk  and  Average  Sample  Number 


! 

S 

f 

For  Truncation  Point  of  2»E (n)  I 

E(n)Q  - mean  of  the  expected  sample  for  each 
test  at  termination  when  input  0*0 q 
using  estimated  parcimeter  values. 

E (n)  - mean  of  the  expected  sample  numbers 

for  each  test  at  termination  when  0=0 
using  estimated  parameter  values. 

Other  column  headings  are  the  seime  as  Appendix  B. 

Tables  are  based  upon  five  hundred  Monte  Carlo 

repetitions. 
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Output  Risk  and  Average  Sample  Number,  Test  One 
=2  Minimum  Sample  Number=5  Asymptotic  0 Statistic 
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[Ut  Risks  and  Average  Sample  Number . Test  Three 
0,=1.5  Minimum  Sample=5  • Exact  G Statistic 
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Appendix  D 
Power 


Tables  are  based  upon  two  hundred  Monte  Carlo 
repetitions  and  are  truncated  at  2»E(n). 


Table  D-I-1 


Power.  Test  One 

0Q=1,  9j^=1.5  Minimum  Sample  Size=5  Asymptotic  0 Statistic 

Power  I k= 


Risks 

CD 

H 

z 

.5 

.75 

1.0 

1.5 

2.0 

3.0 

.20 

.1 

.000 

.000 

.000 

.000 

.000 

.000 

.10 

.000 

.000 

.000 

.000 

.000 

.000 

.20 

.2 

.000 

.000 

.000 

.000 

.000 

.000 

.10 

.000 

.000 

.000 

.000 

.000 

.000 

.20 

.4 

.010 

.000 

.000 

.000 

.000 

.000 

.10 

.015 

.000 

.000 

.000 

.000 

.000 

.20 

.6 

.065 

.045 

.035 

.000 

.000 

.000 

.10 

.035 

.020 

.000 

.000 

.000 

.000 

.20 

.8 

.160 

.135 

.070 

.040 

.010 

.000 

.10 

.105 

.080 

.065 

.010 

.000 

.000 

.20 

1.0 

.350 

.300 

.260 

.180 

.130 

.020 

.10 

.315 

.255 

.235 

.160 

.085 

.045 

.20 

1.2 

.540 

.495 

.330 

.430 

.355 

.425 

.10 

.515 

.500 

.510 

.420 

.465 

.365 

.20 

1.4 

.660 

.695 

.595 

.615 

.595 

.750 

.10 

.795 

.755 

.700 

.700 

.680 

.760 

.20 

1.6 

.845 

.780 

.780 

.820 

.830 

.905 

.10 

.855 

.875 

.835 

.875 

.860 

.900 

.20 

1.8 

.885 

.835 

.870 

.850 

.935 

.965 

.10 

.965 

.880 

.905 

.950 

.935 

.965 

.20 

2.0 

.890 

.900 

.920 

.920 

.955 

1.000 

.10 

.955 

.960 

.945 

.960 

.965 

.995 

.20 

2.2 

.920 

.915 

.925 

.965 

.995 

1.000 

.10 

.995 

.975 

.980 

.980 

.985 

.995 

.20 

2.4 

.940 

.925 

.980 

.980 

.990 

1.000 

.10 

.970 

1.000 

.975 

.990 

.995 

1.000 

.20 

2.6 

.960 

.980 

.940 

.995 

.995 

1.000 

.10 

.980 

.985 

.990 

.985 

1.000 

1.000 

75 


Table  D-I-2 


Risk 

«IN 

.5 

.75 

1.0 

1.5 

2.0 

3.0 

.20 

.1 

.000 

.000 

.000 

.000 

.000 

.000 

.10 

.000 

.000 

.000 

.000 

.000 

.000 

.20 

.2 

.005 

.000 

.000 

.000 

.000 

.000 

.10 

.000 

.000 

.000 

.000 

.000 

.000 

.20 

.4 

.035 

.010 

.000 

.000 

.000 

.000 

.10 

.010 

.005 

.000 

.000 

.000 

.000 

.20 

.6 

.085 

.030 

.005 

.000 

.000 

.000 

.10 

.055 

.040 

.005 

.000 

.000 

.000 

.20 

.8 

.170 

.105 

.030 

.005 

.000 

.000 

.10 

.190 

.075 

.025 

.015 

.000 

.000 

.20 

1.0 

.290 

.205 

.145 

.085 

.025 

.000 

.10 

*.290 

.150 

.150 

.065 

.015 

.000 

.20 

1.2 

.340 

.295 

.250 

.175 

.175 

.050 

.10 

.365 

.280 

.265 

.210 

.080 

.040 

.20 

1.4 

.470 

.470 

.380 

.400 

.345 

.209 

.10 

.460 

.470 

.450 

.420 

.370 

.285 

.20 

1.6 

.525 

.480 

.525 

.505  ■ 

.550 

.615 

.10 

.600 

.590 

.515 

.580 

.585 

.620 

.20 

1.8 

.665 

.635 

.630 

.670 

.760 

.840 

.10 

.690 

.645 

.705 

.730 

.730 

.880 

.20 

2.0 

.725 

.805 

.750 

.775 

.825 

.930 

.10 

.800 

.760 

.795 

.795 

.890 

.945 

.20 

2.2 

.700 

.700 

.750 

.855 

.930 

.965 

.10 

.790 

.805 

.815 

.855 

.955 

.980 

.20 

2.4 

.785 

.790 

.860 

.905 

.930 

.990 

.10 

.860 

.840 

.885 

.915 

.945 

.985 

.20 

2.6 

.785 

.850 

.910 

.955 

.979 

1.000 

.10 

.875 

.865 

.930 

.960 

.965 

.995 

76 


Table  D-II-1 


Power/  Test  Two 

0Q=1/  0j^»1.5  Minimum  Sample  size=5  Asymptotic  G Statistic 

Power  I k= 


Risks 

®IN 

.5 

.75 

1.0 

1.5 

2.0 

3.0 

.20 

.1 

.000 

.000 

.000 

.000 

.000 

.000 

.10 

.000 

.000 

.000 

.000 

.000 

.000 

.20 

.2 

.000 

.000 

.000 

.000 

.000 

.000 

.10 

.000 

.000 

.000 

.000 

.000 

.000 

.20 

.4 

.005 

.005 

.000 

.000 

.000 

.000 

.10 

.000 

.000 

.000 

.000 

.000 

.000 

.20 

.6 

.010 

.010 

.005 

.005 

.000 

.000 

.10 

.005 

.005 

.005 

.000 

.000 

.000 

.20 

.8 

.035 

.045 

.060 

.020 

.000 

.000 

.10 

.015 

.015 

.025  . 

.020 

.000 

.000 

.20 

1.0 

.145 

.155 

.195 

.160 

.095 

.010 

.10 

.045 

.095 

.115 

.115 

.050 

.005 

.20* 

1.2 

.315 

.270 

.390 

.370 

.370 

.240 

.10 

.345 

.250 

.295 

.310 

.300 

.375 

.20 

1.4 

.560 

.490 

.550 

.585 

.615 

.600 

.10 

.655 

.640 

.535 

.645 

.585 

.565 

.20 

1.6 

.620 

.670 

.715 

.715 

.800 

.875 

.10 

.655 

.690 

.725 

.775 

.835 

.895 

.20 

00 

• 

.725 

.760 

.870 

.890 

.810 

.965 

.10 

.800 

.850 

.850 

.905 

.810 

.950 

.20 

2.0 

.880 

.795 

.880 

.900 

.955 

.975 

.10 

.850 

.860 

.810 

.925 

.955 

.985 

.20 

2.2 

.860 

.885 

.850 

.975 

.990 

.990 

.10 

.860 

.945 

.945 

.985 

.980 

1.000 

.20 

• 

CM 

.825 

.895 

.950 

.950 

.990 

1.000 

.10 

.915 

.930 

.955 

.970 

.985 

.995 

.20 

2.6 

.870 

.930 

.965 

.985 

.995 

1.000 

.10 

.900 

.910 

.975 

.990 

1.000 

1.000 

77 


■ 


Table  D-II-2 
Power , Test  Two 


0Q*1,  ®]^*2  Minimum  Sample  Size=5  Asymptotic  G Statistic 

Power  I K= 


Risk 

®IN 

.5 

in 

• 

.20 

.1 

.000 

.000 

.10 

.000 

.000 

.20 

.2 

.000 

.000 

.10 

.000 

.000 

.20 

.4 

.010 

.000 

.10 

.005 

.000 

.20 

.6 

.025 

.101 

.10 

.005 

.000 

.20 

.8 

.080 

.085 

.10 

.040 

.025 

.20 

1.0 

.155 

.220 

.10 

. .090 

.130 

.20 

1.2 

.220 

.335 

.10 

.180 

.190 

.20 

1.4 

.340 

.390 

.10 

.345 

.355 

.20 

1.6 

.435 

.515 

.10 

.435 

.435 

.20 

00 

• 

.540 

.610 

.10 

.580 

.565 

.20 

2.0 

.640 

.655 

.10 

.615 

.680 

.20 

M 

. 

to 

.685 

.735 

.10 

.770 

.815 

.20 

2.4 

.750 

.815 

.10 

.810 

.810 

.20 

2.6 

.771 

.825 

.10 

.840 

.855 

1.0 

1.5 

2.0 

3.0 

000 

.000 

.000 

.000 

000 

.000 

.000 

.000 

000 

.000 

.000 

.000 

000 

.000 

.000 

.000 

000 

.000 

.000 

.000 

000 

.000 

.000 

.000 

005 

.000 

.000 

.000 

000 

.000 

.000 

.000 

055 

.010 

.000 

.000 

025 

.000 

.000 

.000 

150 

.050 

.005 

.000 

115 

.040 

.Cf05 

.000 

280 

.180 

.075 

.000 

170 

.130 

.065 

.000 

300 

.209 

.190 

.070 

305 

.330 

.200 

.085 

485 

.445 

.395 

.260 

425 

.475 

.415 

.280 

690 

.650 

.600 

.570 

675 

.605 

.685 

.535 

685 

.675 

.755 

.830 

715 

.805 

.770 

.770 

770 

.805 

.925 

.895 

775 

.845 

.885 

.875 

850 

.910 

.930 

.980 

855 

.900 

.940 

.965 

875 

.900 

.965 

.985 

875 

.935 

.920 

.975 

78 


Table  D-III-1 


Power / Test  Three 


00=1, 

0j^»1.5 

Minimum  Sample  Size»5 
Power  1 K= 

Exact 

G Statistic 

Risk 

CD 

H 

2i 

.5 

.75 

1.0 

1.5 

2.0 

3.0 

.20 

.1 

.000 

.000 

.000 

.000 

.000 

.000 

.10 

.000 

.000 

.000 

.000- 

.000 

.000 

.20 

.2 

.000 

.000 

.000 

.000 

.000 

.000 

.10 

.000 

.000 

.000 

.000 

.000 

.000 

.20 

.4 

.005 

.005 

.000 

.000 

.000 

.000 

.10 

.000 

.005 

.000 

.000 

.000 

.000 

.20 

.6 

.025 

.030 

.015 

.005 

.000 

.000 

.10 

.010 

.015 

.005 

.000 

.000 

.000 

.20 

.8 

.100 

.100 

.065 

.050 

.010 

.000 

.10 

.045 

.020 

.025 

.035 

.010 

.000 

.20 

1.0 

.315 

.265 

.320 

.125 

.185 

.095 

.10 

.160 

.145 

.175 

.210 

.155 

.110 

.20 

1.2 

.480 

.435 

.535 

.445 

.520 

.480 

.10 

.520 

.525 

.450 

.465 

.515 

.480 

.20 

1.4 

.685 

.625 

.630 

.660 

.705 

.800 

.10 

.765 

.715 

.750 

.720 

.780 

.805 

.20 

1.6 

.765 

.805 

.740 

.820 

.820 

.955 

.10 

.960 

.885 

.870 

.905 

.910 

.930 

.20' 

• 

00 

.885 

.880 

.850 

.920 

.935 

1.00 

.10 

.980 

.965 

.940 

.945 

.955 

.975 

.20 

o 

• 

CM 

.885 

.905 

.925 

.945 

.980 

.995 

.10 

1.000 

.975 

.975 

.965 

.995 

1.000 

.20 

2.2 

.980 

.895 

.955 

.970 

.985 

1.000 

.10 

.990 

.995 

.980 

.995 

1.000 

1.000 

.20 

2.4 

.960 

.955 

.950 

.995 

.990 

1.000 

.10 

.990 

.985 

.990 

.995 

.995 

1.000 

.20 

2.6 

.980 

.965 

.975 

1.000 

.995 

1.000 

.10 

.990 

1.000 

.990 

1.000 

.995 

1.000 

79 


Table  D-III-2 


Power , Test  Three 


00-1. 

0j^»2.O 

Minimxun  Sample  Size=5 
Power  1 K= 

Exact 

G Statistic 

Risk 

®IN 

.5 

.75 

1.0 

1.5 

2.0 

3.0 

.20 

.1 

.000 

.000 

.000 

.000 

.000 

.000 

.10 

.000 

.000 

.000 

.000 

.000 

.000 

.20 

.2 

.000 

.000 

.000 

.000 

.000 

.000 

.10 

.000 

.000 

.000 

.000 

.000 

.000 

.20 

.4 

.005 

.075 

.000 

.000 

.000 

.000 

.10 

.005 

.000 

.000 

.000 

.000 

.000 

.20 

.6 

.060 

.040 

.020 

.000 

.000 

.000 

.10 

.030 

.020 

.010 

.000 

.000 

.000 

.20 

.8 

.120 

.160 

.110 

.030 

.000 

.000 

.10 

.090 

.080 

.060 

.020 

.010 

.000 

.20 

1.0 

.240 

.240 

.210 

.140 

.070 

.025 

.10 

.175 

.195 

.215 

.130 

.050 

.045 

.20 

1.2 

.370 

.375 

.340 

.315 

.340 

.215 

.10 

.265 

.335 

.280 

.345 

.325 

.275 

.20  . 

1.4 

.550 

.530 

.495 

.540 

.510 

.615 

.10 

.475 

.485 

.430 

.520 

.565 

.600 

.20 

1.6 

.610 

.600 

.595 

.690 

.780 

.875 

.10 

.610 

.565 

.660 

.700 

.785 

.860 

.20 

1.8 

.680 

.675 

.750 

.815 

.830 

.955 

.10 

.810 

.730 

.730 

.785 

.890 

.960 

.20 

2.0 

.740 

.725 

.780 

.815 

.915 

.980 

.10 

.850 

.865 

.875 

.880 

.895 

1.000 

.20 

2.2 

.795 

.825 

.815 

.915 

.950 

.990 

.10 

.810 

.900 

.865 

.920 

.980 

.995 

.20 

2.4 

.800 

.800 

.895 

.960 

.990 

1.000 

.10 

.880 

.900 

.945 

.965 

.995 

1.000 

.20 

2.6 

.885 

.845 

.830 

.955 

.990 

1.000 

.10 

.950 

.925 

.940 

.975 

.985 

1.000 

80 


Appendix  E 

AAA 

Comparison  of  Estimates  6,  K,  G 

i 


with  Actual  9 , Ky  G 


I - ■ 

i Table  E-I 


^ /V  A 


Comoarison 

of  Estimates 

9 ULj- 

G with 

Actual  0, 

K,  G 

when  9 = 

1.0 
' ■ A 

A 

Risk 

9 

9 

K 

K 

G 

G 

ASN 

.20 

1.0 

1.16 

.5 

.71 

1.0 

1.11 

17.0 

.15 

1.0 

1.11 

.5 

.70 

1.0 

1.08 

23.2 

.10 

1.0 

1.06 

.5 

.67 

1.0 

1.04 

33.1 

.05 

1.0 

1.03 

.5 

.62 

1.0 

1.02 

52.2 

.20 

1.0 

.98 

.75 

1.06 

1.0 

.98 

9.8 

.15 

1.0 

1.00 

.75 

1.06 

1.0 

1.00 

11.9 

.10 

1.0 

1.03 

.75 

1.03 

1.0 

1.03 

16.7 

.05 

1.0 

.93 

.75 

.96 

1.0 

.93 

22.4 

.20 

1.0 

.95 

1.0 

1.43 

1.0 

.93 

7.4 

.15 

1.0 

.96 

1.0 

1.43 

1.0 

.94 

8.6 

.10 

1.0 

.96 

1.0 

1.36 

1.0 

.95 

11.3 

.05 

1.0 

.94 

1.0 

1.37 

1.0 

.95 

15.0 

.20 

1.0 

.95 

1.5 

2.14 

1.0 

.92 

5.9 

.15 

1.0 

.97 

1.5 

2.09 

1.0 

.95 

6.4 

.10 

1.0 

.95 

1.5 

2.21 

1.0 

.89 

6.8 

.05 

1.0 

.92 

1.5 

2.10 

1.0 

.84 

8.1 

.20 

1.0 

.97 

2.0 

2.93 

1.0 

.91 

5.4 

.15 

1.0 

.96 

2.0 

2.88 

1.0 

.89 

5.6 

.10 

1.0 

.95 

2.0 

2.81 

1.0 

.87 

5.7 

.05 

1.0 

.95 

2.0 

2.77 

1.0 

.87 

6.6 

.20 

' 1.0 

.98 

3.0 

4.29 

1.0 

.92 

5.1 

.15 

1.0 

.99 

3.0 

4.29 

1.0 

.96 

5.1 

.10 

1.0 

.98 

3.0 

4.09 

1.0 

.92 

5.2 

.05 

1.0 

.99 

3.0 

4.22 

1.0 

.96 

5.4 

Table  E-II 


Comparison  of  Estimates  § with  Real 


9 ,k/G  when  6=1.5 


Input 


Risk 

6 

0 

k 

.20 

1.5 

2.12 

.5 

.15 

1.5 

1.96 

.5 

.10 

1.5 

1.99 

.5 

.05 

1.5 

1.86 

.5 

.20 

1.5 

1.79 

.75 

.15 

1.5 

1.86 

.75 

.10 

1.5 

1.80 

.75 

.05 

1.5 

1.78 

.75 

.20 

1.5 

1.69 

1.0 

.15 

1.5 

1.78 

1.0 

.10 

1.5 

1.72 

1.0 

.05 

1.5 

1.68 

1.0 

.20 

1.5 

1.58 

1.5 

.15 

1.5 

1.59 

1.5 

.10 

1.5 

1.58 

1.5 

.05 

1.5 

1.61 

1.5 

.20 

1.5 

1.55 

2.0 

.15 

1.5 

1.53 

2.0 

.10 

1.5 

1.59 

2.0 

.05 

1.5 

1.57 

2.0 

.20 

1.5 

1.49 

2.0 

.15 

1.5 

1.51 

3.0 

.10 

1.5 

1.52 

3.0 

.05 

1.5 

1.53 

3.0 

k 

G 

G 

ASN 

.66 

1.22 

1.64 

46.7 

.64 

1.22 

1.54 

63.4 

.58 

1.22 

1.49 

102.2 

.57 

1.22 

1.42 

145.0 

1.06 

1.35 

1.85 

21.5 

1.01 

1.35 

1.87 

30.1 

1.02 

1.35 

1.82 

42.3 

.95 

1.35 

1.73 

67.1 

1.41 

1.5 

2.10 

15.0 

1.41 

1.5 

2.25 

17.8 

1.36 

1.5 

2.34 

25.6 

1.29 

1.5 

1.95 

37.6 

2.22 

1.84 

2.76 

9.1 

2.12 

1.84 

2.67 

10.3 

2.18 

1.84 

2.71 

14.1 

2.18 

1.84 

2.82 

18.9 

2.89 

2.25 

3.55 

6.7 

2.94 

2.25 

3.49 

8.0 

2.79 

2.25 

3.65 

9.1 

3.03 

2.25 

3.92 

.12.0 

4.52 

3.38 

6.06 

5.7 

4.34 

3.38 

6.55 

6.1 

4.45 

3.38 

6.44 

6.5 

4.56 

3.38 

6.95 

7.3 
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Table  E-III 


Comparison  of  Estimates  9 , K,  G with  Actual  9 , F , G 
when  9=2.0 


Risk 

9 

9 

K 

K 

G 

G 

ASN 

.20 

2.0 

2.99 

.5 

.71 

1.41 

2.18 

17.8 

.15 

2.0 

3.02 

.5 

.70 

1.41 

2.17 

23.8 

.10 

2.0 

2.72 

.5 

.67 

1.41 

1.96 

35.0 

.05 

2.0 

2.73 

.5 

.62 

1.41 

1.86 

52.9 

.20 

2.0 

2.48 

.75 

1.07 

1.68 

2.64 

10.0 

.15 

2.0 

2.40 

.75 

1.06 

1.68 

2.53 

12.6 

.10 

2.0 

2.48 

.75 

1.06 

1.68 

2.62 

17.5 

.05 

2.0 

2.39 

.75 

1.02 

1.68 

2.43 

24.2 

.20 

2.0 

2.21 

1.0 

1.45 

2.00 

3.16 

7.9 

.15 

2.0 

2.24 

1.0 

1.48 

2.00 

3.30 

9.2 

.10 

2.0 

2.17 

1.0 

1.44 

2.00 

3.05 

11.5 

.05 

2.0 

2.24 

1.0 

1.41 

2.00 

3.12 

15.4 

.20 

2.0 

2.09 

1.5 

2.11 

2.83 

4.74 

6.1 

.15 

2.0 

2.06 

1.5 

2.12 

2.83 

4.63 

6.6 

.10 

2.0 

2.11 

1.5 

2.13 

2.83 

4.91 

7.4 

.05 

2.0 

2.11 

1.5 

2.26 

2.83 

5.41 

8.9 

.20 

2.0 

2.02 

2.0 

2.90 

4.0 

7.68 

5.5 

.15 

2.0 

2.01 

2.0 

2.95 

4.0 

7.84 

5.8 

.10 

2.0 

2.04 

2.0 

2.94 

4.0 

8.13 

6.1 

.05 

2.0 

2.09 

2.0 

2.96 

4.0 

8.86 

6.7 

.20 

2.0 

1.98 

3.0 

4.16 

8.0 

17.14 

5.1 

.15 

2.0 

1.96 

3.0 

4.38 

8.0 

19.06 

5.1 

.10 

2.0 

1.99 

3.0 

4.43 

8.0 

21.08 

5.2 

.05 

2.0 

2.00 

3.0 

4.29 

8.0 

19.56 

5.4 

84 


Vita 


Richard  L.  Hoffert  was  born  on  28  December  1940  in 
Caro,  Michigeui.  He  graduated  from  high  school  in  Fremont 
Ohio  in  1958  and  attended  the  United  States  Air  Force 
Academy  from  which  he  received  the  degree  of  Bachelor  of 
Science  in  Public  Policy  in  1962.  Upon  graduation  he  was 
commisioned  as  a 2nd  Lieutenant  in  the  USAF.  He  attended 
pilot  training  and  served  in  a variety  of  flying  and 
logistics  assignments  before  entering  the  Air  Force 
Institute  of  Technology  in  June  1975. 


f 


85 


r 


AD-A035  458 

unclassified 

40A035«e 


AIR  FORCE  INST  OF  TECH  WRIGHT-PATTERSON  AFB  OHIO  SCH— ETC  F/G  12/1 
A monte  CARLO  STUDY  OF  COMPOSITE  SEQUENTIAL  LIKELIHOOD  RATIO  TE— ETC(U) 
DEC  76  R L HOFFERT 


GOR/MA/76-1 


END 

DATE 

FILMED 


NL 


3-77 


SCCumry  classification  of  this  FAGC  r»h«i  Om*  Bnund) 


REPORT  DOCUMENTATION  PAGE 


READ  OfSTRUCnONS 
BEFORE.  COMPLETING  FORM 


2.  GOVT  ACCESSION  NOJ  *■  NECIPIENT’S  CATALOG  NUMGEN 


4.  title  fan^  SuAlltl*> 

A MONTE  CARLO  STUDY  OF  COMPOSITE 
SEQUENTIAL  LIKELIHOOD  RATIO  TESTS  FQJ 
THE  WEIBULL  SCALE  PARAMETER 


S.  TYPE  OF  NEPONT  A PERIOD  COVERED 


S.  PERFORMING  ORG.  REPORT  NUMBER 


T.  AUTHORr«> 

Richard  L.  Hoffert 
Major  USAF 


S.  PERFORMING  organisation  NAME  AND  ADDRESS 

Air  Force  Institute  of  Technology  (AFIT/EN) ^ 
WPAFB  OH  45433  . ^ 


S.  CONTRACT  OR  GRANT  NUMSERT*; 


10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  A WORK  UNIT  NUMBERS 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Air  Force  Flight  Oynanics  Laboratory 
WPAFB  OH  45433 


\ 


12.  REPORT  DATE 

December  1976 


IS.  number  of  PACES 

96 


monitoring  agency  name  a ADORESSfir  MUmnnt  Inm  ContnlUnt  Ottle»)  IS.  SECURITY  CLASS,  (ot  thU  npott) 

Unclassified 


ISa.  OECLASSIFICATION/OOWNCRAOING 
SCHEDULE 


IS.  DISTRIBUTION  STATEMENT  (ot  thio  Koport) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ot  Olo  aAalraet  ontorod  In  Block  20,  II  dittoront  Irom  Hopott) 


e;  lAW  AFR  190-17. 


IS.  KEY  WORDS  (Conthnio  an  rororoo  oldo  II  nocoooory  ond  Idonlllp  by  block  manAarJ 


Weibull  Distribution 
Sequential  Likelihood  Ratio  Test 
Scale  Parameter 
Nuisance  Parameter 

Monte  Carlo  Analysis  ^ 


la  NSUTRAC-T  (Conllnuo  an  ravaraa  oldo  II  nocoocnry  and  Idonlllr  by  block  nuwAarJ  / 

^Three  Sequential  Likelihood  Ratio  Tests  were  constructed  to 
choose  between  two  values  of  the  scale  pareuneter  (6)  of  a Weibull 
distribution  with  a location  parameter  of  zero,  an  unknown  shape 
parameter  (K)  replaced  by  its  maximum  likelihood  estimate,  and 
error  bounds  of  0.05,  0.10,  0.15,  and  0.20. 

The  null  hypothesis  for  all  tests  was  6 equals  1.0.  Two 
alternate  hypotheses  were  used:  6 equals  1.5  and  2.0. 


FORM 
I JAN  71 


EDITION  OF  I NOV  At  IS  OBSOLETE 


UNCLASSIFIED 

security  classification  of  TNIS  page  I 


The  first  two  tests  were  based  upon  a general  procedure 
for  sequential  tests  for  a parameter  in  the  presence  of  a 
nuisance  parameter  by  D.  R.  Cox  in  Sankhya  A,  Vol.  25.  This 
method  replaces  the  likelihood  ratio  with  an  asyit^totically 
equivalent  test  statistic. 

The  tests  were  then  truncated  and  the  effects  of  this 
truncation  upon  the  error  bounds  were  studied tests  wer 
conducted  with  six  input  values  of  K ranging  froBT'rS— to  3.0. 
Newton-Raphson  procedures  were  used  to  estimate  K.  The 
entire  test  was  conducted  using  Monte  Carlo  procedures.  The 
power  of  the  trxincated  tests  was  also  investigated. 

Tabulated  output  data  provides  the  output  Type  One  and 
Type  Two  errors  and  the  average  sample  numbers. 


StCUniTY  CLAUiriCATION  OP  THIS  PAOKfWhM  Oat*  Aitarap) 


